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Multifractal age? 
Multifractal analysis of cardiac interbeat intervals in 

assessing of healthy aging 
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Abstract — 24-hour Holter recordings of 124 healthy people at 
different age are studied. The nocturnal signals of young people 
reveal the presence of the multiplicative structure. This structure 
is significantly weaker in diurnal signals and becomes less evident 
for elderly people. Multifractal analysis allows us to propose 
qualitative and quantitative methods to estimate the advancement 
of the aging process for healthy humans. 

Index Terms — autonomic control, heart rate variability, mul- 
tifractal spectrum 

I. Introduction 

The structure-function-based multifractal analysis relies on 
scaling properties of some statistical measure of the data set 
(see (T)). Namely, if {Xi, i = 1,2, ...,iV} is a time series 
for which the multifractality is investigated, and R(i,n) is a 
function that measures a certain property of a signal in the 
ith box containing n consecutive points of data (boxes do 
not overlap), then the multifractal analysis considers scaling 
properties of all real value moments of R(i,n). Then, the 
question is whether the following power-law dependence on 
box size n exists or not for various real q's: 

Z(n,q)=<\R(i,n)\«> {Xi} ~n T M, (1) 

where < . >{X;} denotes averaging over the corresponding 
data sets. When it is found that such a power-law scaling 
is present for some q, then we say that the studied process 
has fractal structure. Furthermore, if the scaling exponent 
function r(q) is not linear in q then it is said that the process 
is multifractal. The multifractal spectrum of such a process: 
h —> D(h) is obtained from the scaling exponent function 
r(q) by the Legendre transformation (q,r(q)) — > (h,D(h)). 

Different values of q correspond to certain statistical prop- 
erties of a signal. For example, when q = 2, the method 
provides estimates for the Hurst exponent, H = (t(2) + l)/2. 
If q — > oo then we get information about large changes in 
the signal, while q — > — oo characterizes the smooth parts of 
a series. These properties support and motivate our efforts in 
developing methods based on the multifractal approach. 

It often appears that the scaling curves, namely log Z(n, q) 
vs logn plots are, at least at the first sight, linear. Hence, 
according to the above given arguments, the signal can be 
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classified as fractal. However, more careful investigation al- 
lows one to distinguish different scaling regimes. One type of 
scaling is present for the short-time scales and another one 
for the long ones. Such variations of scaling can be related 
to changes in the intrinsic properties of the considered time 
series 0, (6). 

The cardiac interbeat time intervals, so called RR-signals, 
are studied in four generally accepted [7| frequency regimes: 

. high-frequency (HF) (0.15,0.4) Hz, 

. low-frequency (LF) (0.04,0.15) Hz, 

. very-low-frequency (VLF) (0.0033, 0.04) Hz, 

. ultra-low-frequency (ULF) < 0.0033 Hz. 
Such a division is associated with different, physiologically 
distinguishable aspects of cardiac rhythm control. Therefore, 
by considering time scales in the scaling procedure corre- 
sponding to the listed frequency bands, we hope to get insights 
into these aspects. In practical application, assuming that the 
mean cardiac interbeat is 0.8 s, we have assigned size of 
"boxes" for investigation of the scaling properties in signals. 
In particular, we assume n = 32,..., 420 RR-intervals as 
corresponding to the VLF band 

Multifractal analysis based on the study of structure func- 
tions must be applied with care. If the data exhibits noise 
properties then the partial summation must be performed 
before the proper analysis JT). On the other hand, if the data 
has some features of the random walk, there is no need for 
partial summation. Therefore, before multifractal analysis one 
needs to decide whether the data should be treated as a noise 
process or as a random walk. In the following we will show 
what additional information can be gained if the initial data is 
considered both as a noise process and as a random walk. 

The method summarized above was applied to 24-hour 
Holter recordings of cardiac time interbeat signals obtained 
from healthy persons. The obtained results led to the estima- 
tors which allowed us to describe and (at least to some extent) 
quantify the advance of the aging process. The specific details 
are given in Q, Q, |0J. Extensive discussion, numerous 
examples are presented in just cited articles. Here we focus 
only on the essentials. 

II. Methods 

24-hour Holter monitoring of ECG was performed for 124 
healthy subjects: young adults: 36 persons (age 18-26), middle- 
aged adults: 40 persons (age 45-53), elderly: 48 persons (age 
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65-94). Thus, we have analyzed RR-signals consisting of 
time intervals between consecutive normal cardiac contractions 
which originated in the sinus node. 

Multifractal analysis, similarly as in the time-domain statis- 
tics, is found to be robust to data removal (see |9|). However, 
multifractal statistics requires long consecutive series to have 
access to large pool of scales. Additionally, studied series 
should have a stable character. Therefore, we carefully divided 
the obtained signals into two parts corresponding to nocturnal 
and diurnal (daily) activities. The high levels of meanRR, 
SDNN and pNN50 were assumed to specify sleep hours. For 
each RR-signal the period of the six consecutive hours with 
high levels of these quantities was extracted and labeled as 
sleep. Then, the part corresponding to the typical afternoon 
activity was labeled as wake. 

Multifractal analysis relying upon the study of structure 
functions was performed with the well-known approach called 
Wavelet Transform Modulus Maxima (WTMM) [8]. The 
WTMM-based analysis was performed for each person, both 
for wake and sleep parts of his/her signals. This was done 
not only for the proper signal, but also for the correspond- 
ing integrated one. The scaling properties were estimated 
for scales corresponding to the VLF band. Then, after the 
necessary Legendre transformation, the multifractal spectra 
were obtained for wake and sleep parts of every individual 
RR-signal. 




Fig. 1. Assessment of a typical multifractal spectrum: h max - scaling of 
the most dense subset of a signal; Hurst - correlation; hi e f t - extreme, 
and h rig f lt - smooth events scalings; width = /Vight — hi e f t - variety of 
subsets with different scalings 

In Fig. [T] we present the typical multifractal spectrum and 
collect the main parameters characterizing such a spectrum. 
However, real spectra obtained from cardiac time series differ 
considerably from typical ones. Examples of such spectra are 
presented in Fig. (|2]). The evident irregularities (especially in 
the integrated signals) require special care and the correspond- 
ing analysis is quite time-consuming. The details concerning 
such spectra (examples and their discussion) can be found in 

III. Results 

A. A. max in cardiac signals 

The distance A max between h max in the spectrum obtained 
from the integrated signal and h max in a spectrum obtained 
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Fig. 2. Examples of multifractal spectra obtained from RR series. 

in a signal plays a very important role. It allows to distinguish 
additive from multiplicative processes. We have found that 
such a distinction can be made for diurnal and nocturnal parts 
of the signals. The main observation is that the nocturnal 
signals possesses multiplicative character in the VLF region. 
Such organization is quite weak in the diurnal signals. More- 
over, it was found that the multiplicative features are less and 
less pronounced when the age of the investigated individuals 
grows. Hence, weakening of multiplicativity can be considered 
as an indicator of aging. 

B. Multifractal age 

In our analysis we have tested several parameters to char- 
acterize typical multifractal spectra. In particular, for each 
person, we have examined: h max , H (Hurst), hi e f t , h r i g h t , 
A - the width (see Fig{T]for the corresponding illustration) in 
his/hervvoA;e and sleep parts of the signals and their integrated 
counterparts. 

Tedious statistical analysis for sleep and wake parts of the 
signals allowed as to formulate the following five criteria. 

i sleep ^ iwake 

"'max ~ "'max > W 

h wakeMt h large enough(> 1.15), (3) 

£sleep 1 (4) 
^rnax X J vv 

j_j-sleep ^ j_^wake 
H int,sleep is ^ enough (6) 

The notation was explained above and the additional super- 
scripts seem to be self-evident. These criteria were applied to 
all 124 signals obtained from all age groups. 

The results are summarized in the table below. The first row 
tells us how many positive answers were given to the question 
whether each (out of 5) criterion is fulfilled. The main body 
of the table presents the percentage of the number of positive 
answers obtained for three investigated age groups. 





Number of positive answers 


group name: 


1 2 3 4 5 


young 


53 25 11 11 


middle 


3 33 30 20 7 7 


elderly 


10 10 21 29 29 



The number of positive answers clearly increases with age 
and the table has a characteristic skew form. The skewness of 
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the table contents suggests that we can associate the number 
of positive answers to criteria Q - |6]) as indicators of the 
multifractal age. 



IV. Conclusions 

We have applied the multifractal analysis to 124 24-hour RR 
signals divide into wake and sleep parts. The main conclusions 
can be summarized as follows 

• The skewness of the presented table stems from the 
careful multifractal analysis of Holter recordings of RR 
signal obtained from healthy individuals. The number 
of positive answers whether the criteria are met, clearly 
increases with age. Thereby, we suggest that the term 
multifractal age may be useful in the studies of healthy 
aging. 

• The quantity A max leads to distinction between additiv- 
ity and multiplicativity of the investigated RR signals. 
Hence, we suggest that A max should be accepted as an 
additional parameter characterizing multifractality of the 
time series. 
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